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ABSTRACT 

Symmetry adaptation techniques are applied to the determination of the 
intensity of two-photon absorption transitions, between Stark levels of con- 
figurations with opposite parities, for transition ions in finite symmetry en- 
vironments. The equivalence between third order mechanisms and second- 
plus first-order mechanisms is clearly established. A compact formula is 
derived for the intensity of a two-photon transition between two levels with 
well-defined symmetries. A set of (selection) rules which controls the inten- 
sity parameters for inter-configurational transitions is given. The formal- 
ism is illustrated with the case of an ion with Af configuration in tetragonal 
symmetry. 



1. Introduction 

The theory of two-photon absorption spectroscopy for an ion with a partly-filled shell 
embedded in a crystal has received a great deal of attention in recent years [1-13]. It is the 
aim of this paper to develop a formalism for nl — > n'i' inter-configurational two-photon 
transitions when the two involved (one-electron) configurations have opposite parities. This 
problem has been already addressed by Gayen and Hamilton [7], Gayen et al. [7], Leavitt 
[11], Makhanek et al. [4], and Sztucki and Strek [9]. However, little attention has been 
paid to symmetry adaptation techniques (i.e., group-theoretical considerations from both a 
qualitative and a quantitative viewpoint) for inter-configurational two-photon transitions. 
In this paper, we follow the general line adopted in Refs. [12,13] for intra-configurational 
two-photon transitions in order to develop a model incorporating parity violation by the 
odd components of the crystal-field interaction and symmetry adaptation to the chain 
0(3)* D G*. (The group G* denotes the double group of the point symmetry group G of 
the ion site.) The model is applied to the case of the Ce 3+ ion in CaF 2 [7] and LUPO4 
[14]. 

2. Transition element 

We begin with the (electronic) transition matrix element Mj_>/ for a two-photon 
absorption between an initial state i and a final state /. We work in the framework of the 
electric dipolar approximation and use single- mode excitations of energy ftw, wave- vector 
k, and polarization vector S for the radiation field. (For the sake of simplicity, we first 
consider that the two absorbed photons are identical.) Then, the matrix element M^f is 
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given by [15,16] 



Mi 



>f = J2^(f\ 3 -£\ v ) (v\D.£\i) A v = Ei- (E v - hu) (1) 



with evident notations. For linear polarization, the spherical components of £ are 

E = cos 9 £±i = =F— ^= sin6> exp(±z</?) (2) 

v2 

where (6>, <p) are the polar angles of the unit polarization vector £ with respect to the 
crystallographic c-axis. For circular polarization, the components S q (q = — 1,0, +1) of £ 
are given by 

£_ q = -S(q,+l) or £_ q = -S(q,-l) (3) 

for right or left circular polarization, respectively, when the wave vector k is parallel to the 
c-axis. 

The state vectors corresponding to i and / may be (partially) labelled by the irre- 
ducible representation classes of the group G*. We thus have 

\i) = Kzr 7 ) |/) ee \n'£'fv'i) (4) 

where the labels 7 or 7' may serve as multiplicity labels when the dimensions of the 
irreducible representation classes V or V are greater than 1, respectively. The state vectors 
(4) can be expanded in terms of symmetry adapted state vectors as 



\n£ir~f) = \ns£jar~f) c(n£jaT;i) 



(5) 

InY/ry ) = W'sl'j'a'Y'i) c(n'£'fa'V; f) 



with 

\ns£jaT-f) = \ns£jm) (jm\jaT^) 

m 

\n's£'j'a'Y'i) = \n's£'j'm!) (j'm'\j'a'T'i) 



(6) 



where the expansion coefficients of type (JM\JaT^) in (6) are reduction coefficients to 
pass from the {JM} scheme to the {JaF7} scheme [17]. The coefficients c in (5) are 
obtainable from the diagonalization-optimization of the Hamiltonian H describing the ion 
in its environment. (In (5) and (6), s stands for 1/2.) 

The transition matrix element Mi^f given by (1) and (4)- (6) vanishes since the ini- 
tial and final states have opposite parities. In order to generate non-vanishing transition 
matrix elements, we must admit that the parity of the initial and final state vectors is 
not well-defined. This is the case if the crystal-field part of H contains odd components 
either because the site symmetry is non-centrosymmetric or because of the vibration of 
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the ligands. Therefore, we have to replace the initial and final state vectors, respectively, 
by 

\i>=\i)--J2 \n'£'x'T^)(n'£'x'T^\0\n£iT^) (7) 

x' 

and 

1/ > = I/) + ^ Z) HxY'i){nlxY'i\0\n'efY'i) (8) 

X 

In deriving (7) and (8), we have made use of the following points. 

(i) We have supposed that the scrambling between the state vectors of opposite parities 
is restriced to the interaction, via the operator O, of the configurations n£ and n'£'. The 
operator O can be viewed as the odd part of the crystal-field interaction contained in H 
and thus can be developed either as 

k q 

in terms of standard components C*? or as 

k ao 

in terms of G-invariant components C^ ()Tol() . (The label To stands for the identity repre- 
sentation class of the group G.) The one-electron tensor operator C k in (9), with spherical 
components C k (adapted to the chain 0(3) D 0(2)) or symmetry adapted components 
0^ o p o7o (adapted to the chain 0(3) DO), is defined through 

(£'\\C k \\£) = (-Y'[(2£' + l)(2£+l)]^l k Q Q (10) 

(The crystal-field parameters B[kq] in (9a) correspond to the Rajnak-Wybourne param- 
eters in a somewhat different form since we want to reserve the notation B k for the q-th 
component of a tensor B k to be defined below.) 

(ii) Furthermore, we have taken into account that, due to the G-invariance of the 
operator O, the quantum numbers and r'7' remain good quantum numbers for the 
states % and /, respectively. 

(iii) Finally, we have applied the (quasi-closure) approximation (cf. Refs. [18,19]) 

£(n'£Vr 7 ) - E(n£iY<y) = E n , t - E n£ = +A 
E{n£xY'i) - E{n'£'fT'i) = E ni - E n . v = -A ( U ' 

where A = 2nu, a function of the common energy hu of the two identical photons. 
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The next approximation is to suppose that the intermediate states v in (1) are mixed 
by the operator O in a way similar to the one for the states % and / and that, in addition, 
they may have a predominant character of type n£ or n'£' . Consequently, we take 

\v > = \n£x v Y v ^ v ) - — \n'£'x v 'Y v ^ v )(n'£'x v 'Y v ^ v \0\n£x v Y v ^ v ) (12) 

x v 

or 

\v>= \n'£'x v Y v ^ v ) + — \n£x v 'Y v ^ v )(n£x v 'Y v ^ v \0\n'£' x v Y v ^ v ) (13) 

X v 

according to whether as the state v has an n£ or n't' predominant character. 

The transition matrix element (1) now can be non-zero if we replace the initial, inter- 
mediate and final state vectors by the corresponding mixed state vectors. As a result, we 
obtain 

- I A 2 M(2 + l) w = 

{f\0\n£xY'i){n£xY'i\Dj\n'£'x v Y vlv ){n'£'x v Y vlv \D.£\i) 

xx v T 1 v "y v 

+2 (f\D.£\n£x v r vlv )(n£x v r vlv \0\n'£'x v 'r vlv )(n'£'x v 'r vlv \D.£\^ (14) 

x v x v T v "y v 

+ Y (f\D.£\n£x v Y v ^ v )(n£x v Y v ^ v \D.£\n'£'x'Y^)(n'£'x'Y-f\0\i) 

X Xu^v^fv 

where we have neglected the terms in (1/A) 4 in M^f and taken into consideration the 
hermiticity property of the operator O. The notation M(2 + for M^f indicates 

that the transition matrix element has been calculated from second-order time-dependent 
perturbation theory and first-order time-independent perturbation theory. 

It is interesting to note that the result (14) can be obtained in a different way by 
invoking third-order mechanisms and by making use of approximations similar to the ones 
employed for deriving M(2 + l)i— By third-order mechanisms, we mean that we start 
from the transition moment for a three-photon absorption transition where we replace one 
of the operators D. £ by the operator O with change of denominators and appropriate 
permutations. This leads to 

M(3) w = V (flO^iv^D.^iv^DJli) 

— ^12 

V1V2 

— ^12 

V1V2 



(15) 



^12 

V1V2 



The notation M(3)i_>/ for M^f is a reminder that the calculation of the transition matrix 
element arises from third-order time-dependent perturbation theory. The state vectors in 
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(15) are non-mixed state vectors. Furthermore, we assume in (15) that the intermediate 
states v are of either n£ or n'i' type. By using the approximation (11), we get 

Ai 2 (lst sum) = -^A 2 Ai 2 (2nd sum) = - jA 2 Ai 2 (3rd sum) = -^A 2 (16) 

Zi 4r 

so that we finally arrive at the important result 



M(3) W = M(2+1) W (17) 

At this stage, it should be noted that the equivalence described by (17) also holds when the 
two absorbed photons are not identical. In this case, we have to replace A by hio\ + hu)2 
(for more details, see the thesis by one of the authors (M. D.)). 
3. Intensities 

Equation (14) and its analogue in terms of non-identical photons can be worked out 
by using recoupling and symmetry adaptation techniques. In the case of two non-identical 
photons, this leads to 

M i ^ f = (n'£'fT' 1 '\H eff \n£iT 1 ) (18) 
where the effective transition operator H e ff reads 

Heff = £ M[(hk)t] ({{£ 1 £ 2 } h B k } t . U*) (19) 

hkt 

, — » — * 

In eq. (19), {£\£i\ is the tensor product of the polarization vectors £\ and £ 2 for the two 
photons with h = 0, 1, 2. (When £\ = £ 2 , we have only h = 0, 2.) Furthermore, the tensor 
B k in (19) has its spherical components B k given by 

B* = (-)« B[k, -q] (20) 

where the parameters B[kq] are the expansion coefficients of the crystal-field interaction 
O (see eq. (9a)). Hence, k may take the values \t — £'\, \£ — £'\ + 2, ...,£ + £'. Finally, the 
(one-electron) tensor operator U l = u l in (19) is such that 

{t\\v*\\l') = 1 (21) 

and its rank t is restricted by the fact that (£, £', t), (j, j', t), and (/i, k, t) should form three 
triads. 

The parameters M[(hk)t] occurring in (19) can be written as 

M[(hk)t] = ^ [(2t+l)(2h + lp e 2 {n£\r\n'£'f ^WC 1 ^') 2 (^||C fc ||f) 

j h k t\ (1 1 h\ j h k t \ j 1 1 h\ 
{£' £ £){£ £ £'j +Z \£ £' £'){£' £' £) 

(22) 



f h k t ) 

x \ 1 £ £' ) 
I 1 £' £ J 
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where A stands here for Uu\ + Tlojz = E n /£/ — E n i and the factors x, y, and z are 



x = (-Y[i + (-) h ] 



i i 

+ 



!/ = (-)' 



1 +(-)* 1 



z = (-) t+1 y (23) 



It is to be mentioned that, in the special case where the photons 1 and 2 are identical, the 
expressions for M[(hk)t] and H e ff agree with those derived by Leavitt [11] from second 
quantization techniques and by Sztucki and Str§k [9] from ordinary Wigner-Racah calculus. 
(The results of Leavitt and of Sztucki and Strek are derived in a spherical basis rather than 
in a basis adapted to the chain 0(3)* D O*.) 

By making use of symmetry adaptation techniques for the chain 0(3)* DO*, it can 
be shown, from eqs. (5), (18) and (19), that 



= - E E E E E E E E <o'<*' r '; /)* >) 

j'a' ja hkt a" T" 7" a" a 

(2t + I)* (s#|| u *|K7") M[(^fe)t] S[A;a ]* {5i5 2 }^, r „ 7 „ 



/ 



t 



h 



a 1 7 a 1 7 a l o7o 



/ 



3 



t 



3 

ar 7 a'Y'i a"T"-f 



(24) 



where the / coupling coefficients are defined by means of the generic formula [17] 



/ 



h h k 
fJ>i fJ>2 V 



E 

m 1 qrri2 



t\h-mi ( h k j 2 
\—mi q m 2 

0'iTOib'iA»i)* (kq\k/j) (i 2 TO 2 |i2/U2) 



(25) 



with the abbreviations \ii = aiT^i (i = 1, 2) and \i = ar7. (Note that the / coefficient is 
simply a symmetry adapted version of the coupling factor which occurs, besides the reduced 
matrix element, in the ordinary Wigner-Eckart theorem.) In order to obtain eq. (24), we 
have used a property of the / coefficients which arises when /j, = aor 7o (see Ref. [17]). 
The next step is to calculate the intensity strength 



Sr^v = Si ( r)-»/(r') = Yl l M *(r7W(r' 7 ') I 

77' 



(26) 



where the sum on 7 and 7' must be extended over all the Stark components of the initial 
and final states (of symmetry V and T'), respectively. The calculation of Sr^v can be 
achieved by using : (i) eq. (24) twice, (ii) the factorization lemma for the / coefficients 
[17], and (iii) the orthonormality-completeness property for the / coefficients [17]. As a 
net result, we get the compact formula 



= EEE I[hhrfY"-YY'\ {£i£ 2 } r V 7 » ({Si £2}^, 



(27) 



h,h r,r T" 
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The intensity parameters I in eq. (27) may be developed as 
I\hhrfT";TT'} =?«}-' [V] EEEEEE 

j'a' ja j'a' ja, tp tp 

z h (j'aT,jar,hrr'\tpr'')z h Cj'a'r\jarMr''Jpr''r (28) 
£ (^" a ' r ' + ^ r "lja/3r) G'aT' + tpr"|Ja/?r)* 

where the coefficient 

^(jVr'^ar^rr",^") = (2j + l)"* (2*+l)* c(j'a'T';f)*c(jaT;i) {sljWu^sl' j') 

(29) 

turns out to be independent of the label 7". The coefficients ( + | ) in (28) are isoscalar 
factors for the reduction 0(3)* — > 0* [17]. 

The intensity formula (27) has the same form as the one obtained in a previous work 
[12] for intra-configurational two-photon transitions. The intensity parameters I in (27) 
exhibit the symmetry property 

I[hhrfT"; IT'] = I[hhfrT"; FT']* (30) 

which corresponds to the real character of Sr—>r'- Furthermore, they can be factorized 
into two similar factors when the group 0* is multiplicity-free since the sum over (3 in (28) 
disappears in this case. 

The number of a priori independent intensity parameters I in (27) is controlled by 
the property (30) and the following selection rules. 

Rule 1. In most cases, there is no summation in (27) over the multiplicity labels r 
and f. 

Rule 2. The summation in (27) over r" is restricted by the group-theoretical rules 

r" c (h) r" c (h) r"cr® r'* (31) 

where (h) and (h) are irreducible representation classes of 0(3). Similarly, we must have 
r" C (t) and r" C (t) in (28) and (29). 

Rule 3. In the cases where either E\ = £2 or hui = hioii the summation in (27) 
over h and h is restricted to the values 0, 2 but h and h can take the values 0, 1, 2 when, 
simultaneously, S\ 7^ £2 and Tiui 7^ fiU2- 
4. Illustration 

We now illustrate the model developed in the present paper with the case of a 4f — > 5d 
transition for an ion of configuration 4/ 1 in tetragonal symmetry with G = C^ v or D2d- In 
this case, there are four possible transitions since Y and V may be equal to Tq or IV We 
consider only the situation where the two photons are identical. Then, the possible values 
of h and h in (27) are 0,2 and there is no summation on r and f in (27). In addition, the 
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irreducible representation classes V" in (27) are V" = Ai, E for the transitions T 6 — > T 6 or 
r 7 -> T 7 and r" = Si, £ 2 , # for the transitions T 6 -> T 7 or T 7 -> T 6 . By developing (27), 
we get 

S r ^ r =/[02A i; rr] {££}° Ai {{££} 2 Al T + I[2(M i; IT] ({££}° Al )* + 

j[ooa i; rr] |{^}° Al | 2 + j[22A i; rr] |{££} Al | 2 + i[22s ; rr] £ |{££}| y ,| 2 

7" 

(32) 

for T = T' = T 6 or r 7 and 
S T ^ r =/[22S i; rr / ]|{^}|j 2 +/[22S 2 ;rr / ]|{^}|j 2 +/[22E;rr / ]^|{^}| 7 ,,| 2 (33) 



for r 7^ r' = Tq or T 7 . The polarization factors in (32) and (33) can be calculated in an 
0(3) D Coo„ D C4„ D basis to be 

{££}^ = ^i or 

{ ^ i = -L(3cos 2 0-l) or 

1 1 

= ^sin 2 6'cos2^ or + (34) 

{££}| 2 = ^sin 2 0sin2^ or ± i= 



^|{££}| 7 „| 2 = ^sin 2 2# or 



according to as the polarization is linear or circular, respectively. In the detail, equations 
(32) and (33) lead to (for linear polarization) 

Sr 6 ->r 6 = a + bni + cnf + dn 2 
Sr 7 -»r 7 = 0! + b' 7Ti + c' n\ + d! n 2 
Sr 6 -»r 7 = f 7T2 + h 7T 4 + % 7T 5 

r 6 = f 7T2 + h' 7T 4 + l' 7T 5 



(35) 



where the angular functions 7Tj (i = 1, 2, 4, 5) are defined by 

7Ti = 3 cos 2 9 — 1 7T 2 = sin 2 20 7r 4 = sin 4 cos 2 2y> 7r 5 = sin 4 9 sin 2 2</j 

(36) 
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and the various parameters a, • • • , i and a', • • • , i' read 



a = (l/3)/[00A i; r 6 r 6 ] 

b = -(A/3) Re[/[02Ai;r 6 r 6 ] 

c = (l/6)/[22A i; r 6 r 6 ] 

d = (1/2) /[22£ ; r 6 r 6 ] 

/ = (l/2)/[22E;r 6 r 7 ] 

h = (1/2) 7[22s i; r 6 r 7 ] 

i = (l/2)/[22B 2 ;r 6 r 7 ] 



(1/3) /[ooAi ; r 7 r 7 ] 

-(V2/3) Re[/[02Ai;r 7 r 7 ] 

(i/6) i[22A i; r 7 r 7 ] 

(1/2) J[22i?;r 7 r 7 ] 

(1/2) i[22£ ; r 7 r 6 ] 

(1/2) 7[22B i; r 7 r 6 ] 
(1/2) 7[227? 2 ;r 7 r 6 ] 



(37) 



in terms of the intensity parameters 7. The latter parameters are calculable from (28) and 
(29) together with (22). (Equation (28) shows that the parameters 7 can be expressed 
here as the product of two similar factors since there is no internal multiplicity label j3 for 
tetragonal symmetry.) 
4.1. Ce 3 + in CaF 2 

Two- photon experiments for the ion Ce 3+ (of configuration 4/ 1 ) in CaF 2 (with local 
site symmetry G = C^) have been reported in Ref. [7]. The polarization dependence of 
the inter-configurational transition 2 F 5 / 2 (T 7 ) — > 5d(r 7 ) (between the ground state of the 
4/ 1 configuration and the ground state of the 5c? 1 configuration) has been interpreted by 
Gayen and Hamilton [7] on the basis of the group-theoretical formalism from Bader and 
Gold [3] and by Makhanek et al. [4] on a more quantitative basis. 

From equations (35) and (36), the intensity of the 2 F 5 / 2 (r 7 ) — > 5<i(r 7 ) transition can 
be rewritten as 



A + Bsin 2 e + Csm 2 29 



(38) 



with 



A = (1/3) 7[0(L4 i; r 7 r 7 ] + (2/3) 7[22A i; r 7 r 7 ] - 2(^2/3) Re[7[02A i; r 7 r 7 ] 
B = -(1/2) I[22A 1 ; T 7 T 7 ] + V2 Re[7[02A i; T 7 r 7 ]] 

c = -(3/8) i\22A x - r 7 r 7 ] + (1/2) 7[22t; ; r 7 r 7 ]] 



(39) 



The F~ charge compensator in CaF 2 has an equal probability of going into an interstitial 
site along the [100], [010] or [001] directions. Therefore, the intensity strength of the 
2 F 5 / 2 (r 7 ) — >■ 5<i(r 7 ) transition has to be averaged as 



Sr 7 - 



OS([100])rw 7 + S([010])rwr 7 + 5([001])r 7 ^r 7 ) 



Equations (38) and (40) yield 



Sr 7 ^r 7 =A+ (2/3)5 + (2/3)Csin 2 29 



5'r 7 ^r 7 



A+ (1/3)5+ (2/3)C+ (l/3)Bsin 2 6>- (l/3)Csin 2 29 
A+(2/3)B + (2/3)C 



(40) 



(41) 
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for [100] £||[010]), (fc|| [110] ^* || [001]), and (fc||[lll] £||[110]), respectively, which corre- 
spond to the three experimental situations considered by Gayen and Hamilton [7]. The 
intensity formulas (41) formally exhibit the same ^-dependence than those obtained in 
Ref. [7]. Equations (41), where the intensity parameters A, B and C can be a priori cal- 
culated from the formalism developed in section 3, constitutes a justification of the purely 
group-theoretical approach used in Ref. [7]. (Let us remember that the approach used in 
Ref. [7] is based on group-theoretical selection rules for intra-configurational transitions.) 
Finally, it should be noted that the scalar terms (corresponding to h or h = in (27)), de- 
scribed in (41) by the intensity parameters /[OOAi; r 7 r 7 ] and I[02Ai; r 7 r 7 ], do not occur 
in the treatment of Ref. [4] . 
4.2. Ce 3 + in LuP0 4 

Inter-configurational two-photon transitions for Ce 3+ in LUPO4 (site symmetry G = 
D 2 d) have been observed by Piehler [14]. The transitions investigated are of the type 
2 F 5 / 2 (r 6 ) — > 5d (between the ground state of the 4/ 1 configuration and the Stark com- 
ponents of the bd 1 configuration) with linear polarization. Piehler has obtained a first 
transition (at 2hiv = 30460 cm -1 ) whose intensity vanishes when the two (identical) pho- 
tons are polarized along the z-axis and a second one (at 2hu ~ 40000 cm -1 ) whose intensity 
does not vanish for polarization along the z-axis. 

The symmetry species of the bd final states reached in Ref. [14] can be easily deter- 
mined from (35) and (36). Indeed, from (35) and (36), the intensities of the transitions 
2 F 5 / 2 (r6) — > bd(T 6 or T7) can be seen to be 

Sr^re =A' + B' sin 2 9 + C sin 2 29 

Sr 6 ^r r = A" sin 4 9 cos 2 2ip + B" sin 4 9 sin 2 2tp + C" sin 2 29 

where 

A' = (1/3) /[00Ai; r 6 r 6 ] + (2/3) I[22A 1 ; T 6 T 6 ] - 2(^/3) Re[/[02A i; r 6 r 6 ]] 

b' = -(1/2) i[22Ai, r 6 r 6 ] + v^Re[/[02^i; r 6 r 6 ]] 

C = -(3/8) I\22A X - r 6 r 6 ] + (1/2) 7[22£; r 6 r 6 ]] 

a" = (1/2) /[22B i; r 6 r 7 ] b" = (1/2) i[22B 2 - r 6 r 7 ] c" = (1/2) i[22E- r 6 r 7 ] 

Thus, equation (42) shows that only the intensity strength Sr 6 ^r 7 vanishes when the 
electric field is polarized along the z-axis (i.e., 9 = 0). As a consequence, the symmetry of 
the final states is r 7 for the first transition (at 30460 cm -1 ) and Fq for the second transition 
(at 40000 cm -1 ). This result is in accordance with the one obtained by Piehler [14] on the 
basis of selection rules derived from (qualitative) group-theoretical considerations. Here 
again, it is to be mentioned that the intensities (42) are calculable from ah initio principles 
once the intensity parameters (43) are known. 
5. Concluding remarks 

The main result of this paper is formula (27) which gives the polarization dependence 
of an inter-configurational two-photon transition between Stark components arising from 
the configurations n£ and n't' of opposite parities. The obtained formula bears the same 
form as the corresponding one for intra-configurational two-photon transitions. 
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(42) 



The intensity parameters / in (27) are, likewise in the case of intra-configurational two- 
photon transitions, model-dependent. More precisely, they depend on the wave-functions 
for the initial and final states and on the odd crystal-field parameters. They also incorpo- 
rate the information coming from the involved symmetry group. They can be calculated 
in an ab initio way or may be considered as phenomenological parameters. 

Finally, we would like to mention that the model developed in the present work applies 
to the case of n£ N — > n£ N ~ x n' £' two-photon transitions with £ + £' odd. (The passage 
from one-electron configurations to iV-electron configurations does not affect the symmetry 
considerations.) 
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